A Bayesian inference framework to reconstruct transmission trees using epidemiological and genetic data

The accurate identification of the route of transmission taken by an infectious agent through a host population is critical to understanding its epidemiology and informing measures for its control. However, reconstruction of transmission routes during an epidemic is often an underdetermined problem: data about the location and timings of infections can be incomplete, inaccurate, and compatible with a large number of different transmission scenarios. For fast-evolving pathogens like RNA viruses, inference can be strengthened by using genetic data, nowadays easily and affordably generated. However, significant statistical challenges remain to be overcome in the full integration of these different data types if transmission trees are to be reliably estimated. We present here a framework leading to a bayesian inference scheme that combines genetic and epidemiological data, able to reconstruct most likely transmission patterns and infection dates. After testing our approach with simulated data, we apply the method to two UK epidemics of Foot-and-Mouth Disease Virus (FMDV): the 2007 outbreak, and a subset of the large 2001 epidemic. In the first case, we are able to confirm the role of a specific premise as the link between the two phases of the epidemics, while transmissions more densely clustered in space and time remain harder to resolve. When we consider data collected from the 2001 epidemic during a time of national emergency, our inference scheme robustly infers transmission chains, and uncovers the presence of undetected premises, thus providing a useful tool for epidemiological studies in real time. The generation of genetic data is becoming routine in epidemiological investigations, but the development of analytical tools maximizing the value of these data remains a priority. Our method, while applied here in the context of FMDV, is general and with slight modification can be used in any situation where both spatiotemporal and genetic data are available

Evolution of foot-and-mouth disease virus intra-sample sequence diversity during serial transmission in bovine hosts

RNA virus populations within samples are highly heterogeneous, containing a large number of minority sequence variants which can potentially be transmitted to other susceptible hosts. Consequently, consensus genome sequences provide an incomplete picture of the within- and between-host viral evolutionary dynamics during transmission. Foot-and-mouth disease virus (FMDV) is an RNA virus that can spread from primary sites of replication, via the systemic circulation, to found distinct sites of local infection at epithelial surfaces. Viral evolution in these different tissues occurs independently, each of them potentially providing a source of virus to seed subsequent transmission events. This study employed the Illumina Genome Analyzer platform to sequence 18 FMDV samples collected from a chain of sequentially infected cattle. These data generated snap-shots of the evolving viral population structures within different animals and tissues. Analyses of the mutation spectra revealed polymorphisms at frequencies >0.5% at between 21 and 146 sites across the genome for these samples, while 13 sites acquired mutations in excess of consensus frequency (50%). Analysis of polymorphism frequency revealed that a number of minority variants were transmitted during host-to-host infection events, while the size of the intra-host founder populations appeared to be smaller. These data indicate that viral population complexity is influenced by small intra-host bottlenecks and relatively large inter-host bottlenecks. The dynamics of minority variants are consistent with the actions of genetic drift rather than strong selection. These results provide novel insights into the evolution of FMDV that can be applied to reconstruct both intra- and inter-host transmission routes

Transmission dynamics and prospects for the elimination of canine rabies

Rabies has been eliminated from domestic dog populations in Western Europe and North America, but continues to kill many thousands of people throughout Africa and Asia every year. A quantitative understanding of transmission dynamics in domestic dog populations provides critical information to assess whether global elimination of canine rabies is possible. We report extensive observations of individual rabid animals in Tanzania and generate a uniquely detailed analysis of transmission biology, which explains important epidemiological features, including the level of variation in epidemic trajectories. We found that the basic reproductive number for rabies, R<sub>0</sub>, is very low in our study area in rural Africa (∼1.2) and throughout its historic global range (<2). This finding provides strong support for the feasibility of controlling endemic canine rabies by vaccination, even near wildlife areas with large wild carnivore populations. However, we show that rapid turnover of domestic dog populations has been a major obstacle to successful control in developing countries, thus regular pulse vaccinations will be required to maintain population-level immunity between campaigns. Nonetheless our analyses suggest that with sustained, international commitment, global elimination of rabies from domestic dog populations, the most dangerous vector to humans, is a realistic goal

One Health contributions towards more effective and equitable approaches to health in low- and middle-income countries

This research was supported by the UK Biotechnology and Biological Sciences Research Council (BB/J010367/1) and the UK Zoonoses and Emerging Livestock Systems Initiative (BB/L017679/1, BB/L018926/1 and BB/L018845/1) (S.C., J.E.B.H., J.S., J.B., A.D., J.A.C., W.A.d.G., R.R.K., T.K., D.T.H., B.T.M., E.S.S., L.W.). The Wellcome Trust provided supported for K.H. and A.L. (095787/Z/11/Z) and K.J.A. (096400/Z/11/Z). The US National Institutes of Health provided support for J.A.C. (R01AI121378) and M.P.R. (R01AI121378, K23AI116869).Emerging zoonoses with pandemic potential are a stated priority for the global health security agenda, but endemic zoonoses also have a major societal impact in low-resource settings. Although many endemic zoonoses can be treated, timely diagnosis and appropriate clinical management of human cases is often challenging. Preventive ‘One Health’ interventions, e.g. interventions in animal populations that generate human health benefits, may provide a useful approach to overcoming some of these challenges. Effective strategies, such as animal vaccination, already exist for the prevention, control and elimination of many endemic zoonoses, including rabies, and several livestock zoonoses (e.g. brucellosis, leptospirosis, Q fever) that are important causes of human febrile illness and livestock productivity losses in low- and middle-income countries. We make the case that, for these diseases, One Health interventions have the potential to be more effective and generate more equitable benefits for human health and livelihoods, particularly in rural areas, than approaches that rely exclusively on treatment of human cases. We hypothesize that applying One Health interventions to tackle these health challenges will help to build trust, community engagement and cross-sectoral collaboration, which will in turn strengthen the capacity of fragile health systems to respond to the threat of emerging zoonoses and other future health challenges. One Health interventions thus have the potential to align the ongoing needs of disadvantaged communities with the concerns of the broader global community, providing a pragmatic and equitable approach to meeting the global goals for sustainable development and supporting the global health security agenda.Publisher PDFPeer reviewe

Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)

Relationship of edge localized mode burst times with divertor flux loop signal phase in JET

A phase relationship is identified between sequential edge localized modes (ELMs) occurrence times in a set of H-mode tokamak plasmas to the voltage measured in full flux azimuthal loops in the divertor region. We focus on plasmas in the Joint European Torus where a steady H-mode is sustained over several seconds, during which ELMs are observed in the Be II emission at the divertor. The ELMs analysed arise from intrinsic ELMing, in that there is no deliberate intent to control the ELMing process by external means. We use ELM timings derived from the Be II signal to perform direct time domain analysis of the full flux loop VLD2 and VLD3 signals, which provide a high cadence global measurement proportional to the voltage induced by changes in poloidal magnetic flux. Specifically, we examine how the time interval between pairs of successive ELMs is linked to the time-evolving phase of the full flux loop signals. Each ELM produces a clear early pulse in the full flux loop signals, whose peak time is used to condition our analysis. The arrival time of the following ELM, relative to this pulse, is found to fall into one of two categories: (i) prompt ELMs, which are directly paced by the initial response seen in the flux loop signals; and (ii) all other ELMs, which occur after the initial response of the full flux loop signals has decayed in amplitude. The times at which ELMs in category (ii) occur, relative to the first ELM of the pair, are clustered at times when the instantaneous phase of the full flux loop signal is close to its value at the time of the first ELM

Cross-disciplinary demands of multihost pathogens

The dynamics of infectious disease spread depend on host population contact structure. Heterogeneities in this contact structure can arise from various forms of demographic and spatial phenomena. Craft et al. (this issue) have constructed an exploratory simulation model of the spread of canine distemper virus through a multispecies carnivore community. Each species in this community is modelled with a contact structure reflecting host social organization, ranging behaviour, and likely interspecific contact patterns. The results are used to infer the possible roles of different species in determining the observed spatio-temporal incidence of canine distemper virus in Serengeti lions during an outbreak in 1993-9

Reconstructing the origin and transmission dynamics of the 1967–68 foot-and-mouth disease epidemic in the United Kingdom

A large epidemic of foot-and-mouth disease (FMD) occurred in the United Kingdom (UK) over a seven month period in Northwest England from late 1967 to the summer of 1968. This was preceded by a number of smaller FMD outbreaks in the country, two in 1967, in Hampshire and Warwickshire and one in Northumberland during 1966. The causative agent of all four events was identified as FMD virus (FMDV) serotype O and the source of the large epidemic was attributed to infected bone marrow in lamb products imported from Argentina. However, the diagnostic tools available at the time were unable to entirely rule out connections with the earlier UK FMD outbreaks, as well as other potential sources from Europe. The aim of this study was to apply molecular sequencing to investigate the likely source of this epidemic using VP1 region and full genome (FG) sequences determined directly from clinical epithelium samples (n = 13) or cell culture isolates (n = 6), from this and contemporary outbreaks in the UK, Europe and South America. Analysis of the VP1 sequences provided evidence for at least three separate incursions of FMDV into the UK including one independent introduction that was responsible for the main 1967/68 epidemic. Analysis of FG sequences from the main 1967/68 outbreak (n = 10) revealed nucleotide substitutions at 94 genomic sites providing evidence for the linear accumulation of nucleotide substitutions (rate = 2.42 × 10−5 nt substitutions/site/day). However, there were five samples where this linear relationship was absent, indicating evolutional dormancy of the virus, presumably outside a host. These results help define the evolutionary dynamics of FMDV during an epidemic and contribute to the knowledge and understanding from which to base future outbreak control strategies